Bose-Einstein condensation in multilayers 
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The critical BEC temperature Tc of a non interacting boson gas in a layered structure like those of 
cuprate superconductors is shown to have a minimum Tc,m, at a characteristic separation between 
planes am- It is shown that for a < Om, Tc increases monotonically back up to the ideal Bose gas 
To suggesting that a reduction in the separation between planes, as happens when one increases the 
pressure in a cuprate, leads to an increase in the critical temperature. For finite plane separation 
and penetrability the specific heat as a function of temperature shows two novel crests connected 
by a ridge in addition to the well-known BEC peak at Tc associated with the 3D behavior of the 
gas. For completely impenetrable planes the model reduces to many disconnected infinite slabs for 
which just one hump survives becoming a peak only when the slab widths are infinite. 
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Since London first suggested^ that superfluidity in liq- 
uid '^He might well be a manifestation of Bose-Einstein 
condensation (BEC) of the helium atoms before in- 
teratomic interactions are "svifitched on," BEC in lay- 
ered systems began to be studied to understand helium 
films^i^i^. The discovery of high- Tc superconductivity 
stimulated renewed interest in compounds with layered 
structures^ in which a BEC mechanism seems to be an 
essential feature to explain high critical temperatures®. 
The so-called "Uemura plot"— of data from muon-spin 
relaxation (/iSR), neutron and Raman scattering, and 
angle-resolved photoemission (ARPES) measurements 
exhibits Tc vs Fermi temperatures Tp = Ep/ks where 
Ep the Fermi energy and fc^ the Boltzmann constant. 
Empirical TcS of many cuprates straddle a line parallel to 
the Uemura-plot diagonal line associated with the simple 
BEC formula Tq ~ SMffri^j^^ /mkB ^ Q.2l8Tp corre- 
sponding to an ideal gas of bosons of mass m ~ 2m* and 
number density ub — na/2 where m* is the individual- 
charge-carrier effective mass and Us their number den- 
sity. The parallel line of data is shifted down from Tq by 
a factor of 4-5. This has been judged^^ a "fundamental 
importance of the BEC concept in cuprates." In addi- 
tion, the possibility of creating BECs or superfluidity of 
ultracold fermions^ in optical latticesi^, along with the 
expected observation of BEC of excitons (electron-hole 
pairs) in semiconductorsii, have further revived theoret- 
ical and experimental^^'^'^ efforts to better understand 
the behavior of quantum gases in layered geometries. 

Most models based on layered structure o®'^^ simulat- 
ing quasi-2D high- Tc superconductors, or other to study 
BEC'^iiiii^, rely on a single-boson hopping interaction 
term producing nearest-interlayer couplings in one spa- 
tial dimension while moving freely in the other two direc- 
tions. The energy spectrum is typically of the form ek = 
h^{kl + kl)/2m + ek, with efc^ = {h'^/Ma'^){l - cosfc^a) 
where a is the plane separation and the constant h? / Ma'^ 
is a measure of the bosonic Cooper pair hopping prob- 
ability between planes. In the case of Cu02 planes in 
cuprate superconductors a boson effective mass M in the 



z-direction is used with m the mass in the x-, y-directions 
while M > m reflects cuprate anisotropy. For fc^a <C 1 
one recovers the expected result ek, — lT?k^/2M . 

Motivated by possible applications to quasi-2D super- 
conductors or to ^He films, in this Brief Report we discuss 
critical temperature Tc and the specific heat Cy results 
in a non-interacting boson gas of N particles within an 
infinite multilayered structure of equally-spaced planes of 
variable penetrability stacked in the ^-direction. 

The system of parallel planes is simulated with an ex- 
ternal periodic delta potential of the ID Kronig-Penney 
(KP) type along the z-direction while bosons are free in 
the two remaining directions. Two important differences 
with previous studies^ '^^i"'^'^'^^ are: a) the bosons are al- 
lowed to move over all space through permeable planes 
instead of constraining them to move only over the plane 
surfaces, and b) the anisotropy expected from a layered 
material is introduced naturally by the delta KP poten- 
tial. In addition, the delta strengths avoid the artifi- 
cial introduction of large masses and hopping parame- 
ters in the z-direction to tune the mobility or tunneling 
across planes. At very low energies the KP potential re- 
duces to the energy expression ek, used in Refsi^ i^^'^^i^^ 
which only takes into account the lowest allowed energy 
band. If Tq is the critical BEC temperature of the free 
3D boson gas and the associated thermal wavelength 
is Aq = h / \/2TrmkBTo a diniensionless "impenetrabil- 
ity" Pq = PXq ^ of layers is introduced in terms of 
the KP penetrability parameter P [see ((TJ below]. Thus, 
Pq = implies perfect layer transparency and Pq — oo 
fully opaque layers. In the latter limit, our model re- 
duces to that of an infinite number of uncoupled slabs of 
thickness a and infinite lateral extent. The specific heat 
has been obtained in these systems to model *He thin- 
film propertiesi2i^iS^i2^ and finite-size effects on BEC^''-. 
In these infinite slabs the BEC critical temperature is 
zero but the specific heat shows a smooth maximum at 
a temperature that depends on the slab thickness. Some 
authors22i2^i2i have associated this "hump" with a BEC 
signature. However, we show here that for finite Pq this 
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hump coexists with a nonzero critical BEC temperature 
in addition to a second maximum and when Pq oo the 
BEC temperature goes to zero whereas the hump per- 
sists to become the famihar BEC with a cusped sharp 
peak only in the limit a/Ao oo. 

The Schrodinger equation for a boson of mass m 
is separable in the x, y and z-directions so that the 
single-particle energy as a function of the wavevector 
k = {kx,ky,k^) is simply Sk = Sk^ + Sk^ + e^^, where 
£k^,ky = h^kl y/2m with kx,y = 2'Knx,y/L and nx,y 
0, ±1, ±2, • • • , i.e., in the x and y-directions particles arc 
free and we assume periodic boundary conditions in a 
box of size L. In the z-direction the particles are subject 
to the KP periodic delta (or "Dirac comb") potential in 
ID'^i and the energies £k^ are implicit in the transcen- 
dental equation 
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where a^ = 2mek^/h^ and fi^P/m is the delta-interaction 
strength or equivalently, the layer impenetrability pa- 
rameter. The energy £fe^,j depends on the parameters P 
and a where the band structure is labeled with the index 
j = 1, 2, ... The trivial free-particle energy dispersion re- 
lation in the z-direction, Sk^j h^k1/2m, is recovered 
in the limit P — > 0, while P oo yields sin(Q;a) — > 
which corresponds to a system of confined bosons inside 
a semi-infinite slab of width a, a situation extensively 
discussed in the literature (seeSi and references therein). 
For aa <C 1, (HI can be expanded as 



Pa + (1 - cos k^a) 
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where the last expression is the dispersion relation used 
^q6 ,16^18 _ Pqj. gjjj^g bosons tunnel through the layers 
as a more realistic model would demand. 

The grand potential il{T, V, fi) = —pV (with p the 
pressure, V the volume and fi the chemical potential) 
leads to all the system thermodynamic properties. In our 
multilayered system it becomes, after integrating over kx 
and fc,. 
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where z = e*^^ is the fugacity, /3 = l/fc^T, and we have 
explicitly separated the contribution JIq = ksT ln[l — 
ze~^'^°] of the lowest-energy state Eq. Here ga{t) stands 
for the Bose function^^. The integration is over the first 
Brillouin zone and the sum is over the allowed bands. 
Note that the ground-state energy depends on both P 
and a in ([1]). 

The average number of particles N — — [dil/ dfi)rp y 
for a free boson gas in the multilayered configuration is 



then 
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where No{T) = [z ^e^^" — 1] "'^ is the number of parti- 
cles condensed in the ground state Eq while the last term 
in ^ is the number of particles in excited states. As 
usual, the BEC critical temperature Tc is obtained by 
finding the temperature below which the fractional num- 
ber of particles in the ground state just ceases to be neg- 
ligible upon cooling and the chemical potential reaches 
its maximum value Eq > making the fugacity equal 
zq = 6^"^° . For convenience and clarity we use the same 
particle density N/V in the layered system as that for 
the free ideal Bose gas (IBG) with the BEC critical tem- 
perature To- Thus fc^To serves as an energy scale while 
Ao = h / \/2jrmkBTo provides a length scale. 

The effects of the KP layers on BEC are shown in Fig. 
[T] where Tc/Tq decreases as Pq increases for different val- 
ues of a. Note that the critical temperature is very sen- 
sitive to the number of energy bands intervening in the 
numerical calculations. Here, it is sufhcient to include 
10 bands to achieve convergence for the values of a and 
Pq shown in the figure. When Pq = (perfect barrier 
transparency) we recover the IBG results. In the oppo- 
site limit Po — > oo (impenetrable barriers) Tc/Tq vanishes 
although strictly speaking the model does not become a 
true two-dimensional system. 
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FIG. 1: Critical temperature Tc in units of To as a function 
of Po for different values of o/Aq. 

A novel result shown in Fig. [2] is that for finite Pq, 
Tc/Tq diminishes from 1 down to a minimum value as 
the plane separation a decreases. For Pq spanning at 
least five orders of magnitude (see Table |T| the critical- 
temperature minima show a linear dependence with a^n 
so that Tcjn = (To/Ao)a„i or (Ao/Ac,m)^ = Om/Ao, where 
the subscript m means minimum and Xc.m is the thermal 
wave length associated to Tc.m- Since — > and 



TABLE I: Values of plane separation Um for which the critical 
temperature reaches its minimum for different Pq. 



Po 


10 


10^ 


10^ 


10^ 


10^ 


am/ \) 
Tcm 1 To 


0.874 
0.927 


0.550 
0.572 


0.314 
0.337 


0.177 
0.191 


0.100 
0.107 



flm — >■ for Pq ~* oo this linear relation seems to hold 
for all Po > 10. For Po< 10, T^,m/To approaches 1 while 
am oo. Surprisingly, further reduction in a brings 
an increase in Tc/Tq which asymptotically reaches unity. 
This can be understood from the KP dispersion relation 
^ since for small values of a sin(Q;a)/aa ~ 1 and Sk^j — 
lT?kl/2m so that one recovers the 3D IBG regime when 
a ^ 0, as expected. We also note in Fig. [2] that Tc/Tq 
is not symmetric wrt the minimum value. For a/Ao ^ 1, 
Tc/Tq increases slower than for a/Ao ^ 1 which reflects 
the plane influence on the bosons. We also note that for 
a/Ao 3> 1 one must include more than ten energy bands 
to achieve convergence in the numerical calculations. 
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FIG. 2: Critical temperature as a function of a/Ao for different 
values of Pq. For a S> Ao the critical temperature is very 
sensitive to the number of energy bands used in the numerical 
calculations. Thin, dashed and normal lines correspond to 
10, 100 and 1000 bands, respectively. To attain convergence 
we need to consider more than ten bands (thin line) which 
contrasts with the behavior in the a < Ao where the inclusion 
of few bands is enough to obtain convergence 



The specific heat C'v = -T [d^n/dT^]^^ is 
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where 
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For T < Tc the chemical potential /i — Eq is a, con- 
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FIG. 3: 3D plot of the specific heat as function of temperature 
and plane separation a for Po = 100. 
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FIG. 4: Specific heat per particle as a function of temperature 
for several Po and a/Ao = 0.1. The two-dimensional character 
of the layered structure is easily seen as Po increases. The 
arrow indicates the initiation of the IBG regime 

stant so that dfi/dT = which simplifies the equation 
for Cv/NkB- 

Numerical results are summarized in Fig. [3] of 
Cv/NkB vs. logT/To and vs. log a/Ao for an interme- 
diate dimensionless penetrability Po = 100. In contrast 
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with the IBG case, in addition to the famihar BEG peak 
for finite Pq two humps appear that are connected by a 
ridge when T > Tc] ridge heights depend on the plane 
separation. When a > a„i only one hump survives and 
becomes the EEC cusped peak in the limit a — > oo inde- 
pendently of Pq. Our model shows the distinct nature of 
the ever present BEG critical temperature plus a second 
characteristic temperature where Cy/NkB is maximum. 
In previous models (here equivalent to Pq oo this 
latter temperature is identified as a signature of the BEG 
critical temperature. However, our results indicate that 
this characteristic temperature should be associated with 
a different feature, namely an accumulation of bosons^^ 
in the lowest energy levels. 

While Cv/Nks of the IBG decreases monotonically 
for T > Tc = To, in the present model there is a re- 
gion defined by a < am where Cy/NkB varies non- 
monotonically with T having a set of local minima (the 
ridge in Fig. [3]) where A ~ 2a. In this case, the single- 
particle (thermal) average wavefunction nodes coincide 
with the plane locations resulting in the particle motion 
freezing out in the ^-direction. In particular, for a — 
the minimum of Cv{T) / Nks attains its lowest value ~ 1 
so that the "ridge" reveals a closer resemblance to a 2D 
system. As Pq increases this 2D-like feature dominates 
as the ridge levels out over a broader temperature region 
(see Fig. IHfor a/Ao = 0.1). In addition, for temperatures 
below Tc the value A ~ 2a marks a functional crossoever 
of specific heat from the standard 3D IBG Cy cx T^f^ to 
the 2D linear Cy oc T behavior. For sufficiently high tem- 
peratures, layer effects are negligible regardless of separa- 



tion a or of penetrability Pq. In this regime Cv{T) /Nks 
approaches the IBG value. Since Cy attains this value 
at the same T (implying the same A) we can identify an 
a-dependent correlation length that manifests itself when 
A ~ 0.8a as marked by the arrow in Fig. [4]for a/Ao = 0.1 
and T/To ~ 170i^. 

To conclude, we have shown that the BEG-like critical 
temperature as well as the specific heat of a boson gas be- 
tween penetrable and periodically-spaced planes change 
drastically compared to previous models based on hop- 
ping, Hubbard^ or slab^^ configurations, or when there 
is no confining geometrical structure whatsoever. Even 
though analytical approximations for the boson disper- 
sion relation at low energies commonly used to describe 
layered systems have proved to be a useful guide to give 
more accurate results, they are significantly improved 
when the more accurate dispersion relation including 
more than one energy band is used. The appearance of 
two characteristic temperatures associated with the max- 
ima of Cv {T) in addition to the BEG transition tempera- 
ture due to the confinement parameters may explain the 
variation of Tc with pressure in high-Tc superconductors 
the transition is related with a BEG-like phenomenon. In 
the limit Pq oo (opaque or perfectly-decoupled infinite 
slabs of width a) the BEG-like transition temperature 
vanishes leaving only one maximum which becomes the 
usual cusped-peak BEG transition only when a 3> Aq. 
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